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ABSTRACT 

We analyze the fundamental plane projections of elliptical galaxies as a function of lumi- 
nosity, using a sample of « 80, 000 galaxies drawn from Data Release 4 (DR4) of the Sloan 
Digital Sky Survey (SDSS). We separate brightest cluster galaxies (BCGs) from our main 
sample and reanalyze their photometry due to a problem with the default pipeline sky subtrac- 
tion for BCGs. The observables we consider are effective radius (R e ), velocity dispersion (a), 
dynamical mass (Md yn oc R e cr 2 ), effective density (<r 2 /R^), and effective surface brightness 
(fj, e )- With the exception of the L-Md yn correlation, we find evidence of variations in the slope 
(i.e. the power-law index) of the fundamental plane projections with luminosity for our normal 
elliptical galaxy population. In particular, the radius-luminosity and Faber-Jackson relations 
are steeper at high luminosity relative to low luminosity, and the more luminous ellipticals 
become progressively less dense and have lower surface brightnesses than lower luminosity 
ellipticals. These variations can be understood as arising from differing formation histories, 
with more luminous galaxies having less dissipation. Data from the literature and our reanaly- 
sis of BCGs show that BCGs have radius-luminosity and Faber-Jackson relations steeper than 
the brightest non-BCG ellipticals in our sample, consistent with significant growth of BCGs 
via dissipationless mergers. The variations in slope we find in the Faber-Jackson relation of 
non-BCGs are qualitatively similar to that reported in the black hole mass-velocity dispersion 
(Mbh-c) correlation. This similarity is consistent with a roughly constant value of Msa/M* 
over a wide range of early type galaxies, where M* is the stellar mass. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: fundamental parameters - galaxies: statistics 



1 INTRODUCTION 

Early-type galaxies are observed to populate a tight plane - the 
fundamental plane (FP) - in the space defined by their effective 
radii, velocity dispersions, and surface brightnesse s (Dressier et al. 
19871 : iDiorgovski & Davisll987l : lBender et alJl992l ; lBurstein et all 
1997h . Such a correlation is expected theoretically if early-type 
galaxies are in virial equilibrium. The observed plane is, how- 
ever, tilted with respect to the simplest virial theorem expectation, 
which must be accounted for in galaxy formation models. This tilt 
could be caused by a variation in the dynamical mass-to-light ra- 
tio for ellip tical galaxies as a result of a varying dark matter frac- 
tion (e.g.. IPadmanabhan et al]|2004l: iBovlan-Kolchin et alj|200ll) 
or stellar population variations (e.g. jGerhard et alJuOOlf) . Another 
explanation for the FP tilt could be non-homology in the surface 
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brightness profiles of ell iptical galaxies (e.g., iGraham & Collessl 
ll997l : lTruiilloetal.ll2004h . 

Projections of the FP contain additional information about the 
properties of early typ e galaxies. These includ e, e.g., the Faber- 
Jackson (FJ) relation dFaber & JacksorJ Il976[) between velocity 
dispersion and luminosity and the radius-luminosity relation. Al- 
though the FP itself is quite tight ov er a wide range of early type 
galaxies (e.g.. iBernardi et aT, 2003a), there are a number of indi- 
cations of variations in the FP projections themselves with galaxy 
luminosity (which we refer to as "variations in slope" or "curva- 
ture" in the FP projections). For example, the slope of the Faber- 
Jackson relation decreases in lower luminosity ellip ticals jTonrvl 
| l98ll:lDavies et alJI 19831: iMatkovic & Guzmanll200a) . In addition, 
lOegerle & Hoessel J 1 99 lh find that while brightest cluster galaxies 
(BCGs) lie on the same FP as lower mass ellipticals, their central 
velocity dispersions are roughly constant with increasing luminos- 
ity and their effective radii increase steepl y with luminosity, si gnif- 
icantly more so than f or normal ellipticals. lLauer et alj d2006ah and 
Bernardi ct al. (2006) have confirmed these differences in the FP 
projections of BCGs with larger data sets. These trends can be un- 
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derstood if dissipationless (gas-po or) mergers become increasi ngly 
important for luminous ellipticals dBovlan-Kolchin et ai]|2006h . 

In addition to global correlations among the galaxy prop- 
erties themselves, the masses of central black holes in early- 
type galaxies correlate well with the veloc i ty dispersions 
dFerrarese & Merrittl I200M iGebhardt et alj I2OO0I; iFerrarese et~aT 



2001; Tre maine et alj|2002f). bulge concentratio ns ( Graham et al 
2001), and stellar masses dMagorrian et al.l l 1998; Marconi & Hun 
Flaring & Rix 2004) of their host galaxies. IWvitb 



2003; 



and lGreene & Hoi ( 120060 have argued that the Mw — a relation 
steepens at high black hole masses and flattens at low masses, re- 
spectively. This is similar to the luminosity dependence in the FJ 
relation noted above. If supported by larger data sets, the similar 
luminosity dependence of the A/bh — o and FJ relations could have 
important implications for understanding the relationship between 
the formation of galaxies and their central massive black holes. 

The purpose of the present paper is to quantify the variation 
in the FP projections with luminosity using the exquisite statis- 
tics made possib le by the Sloan Digital Sky Su rvey (SDSS) Data 
Release 4 (DR4; lAdelman-McCarthv et al.l l2006). These variations 
have, until now, been largely overlooked by previous FP studies, be- 
cause simple power-law assumptions were often made when fitting 
the data. This can be seen, for insta nce, in the sample of ~ 9000 
galaxies from lBernardi et al] d2003bl) . where hints of the variations 
we discuss are quite noticeable in the plots. In the next section we 
describe the data used in our analysis (Sj2j. We also describe a prob- 
lem in the standard pipeline SDSS photometry of BCGs that leads 
us to reanalyze the BCGs separately from the bulk of the "nor- 
mal" elliptical galaxy population. Having defined our sample, in 
ij3] we examine the observed FP projections and quantify curva- 
ture in these projections as a function of luminosity. We also use 
published photometry of nearby (non-SDSS) BCGs to discuss the 
relationship between our results and the properties of BCGs. In ij4] 
we discuss our results and their connection to galaxy-BH correla- 
tions. Throughout this paper, we use a standard cosmology with 
(Q M , Oa, h) = (0.3,0.7,0.7). 



2 DATA 

The data used in this paper wer e obtained from the NYU Value 
Added Galaxy Catalog (VAGC) (Blanton et al.ll2005h whic h uses 
the SDSS Data Release 4 dAdelman-McCarthv et al]|2006t) . This 
includes both the im aging and spectroscopic catalogs, as well as the 
ISersiddl963lll968h model fitting parameters dBlanton et al.l2003bh 
in 5 ban ds (ugriz). We also us e the VAGC if-corrections (version 
4_1_4) dBlanton et al . 2003a). SDSS magnitudes are calibrated to 
the AB system dOke & Gunnlll983h . in which a magnitude object 
has the same counts as a F„ = 3631 Jy source (this zero point 
has been confirmed in the r band). We use SDSS magnitudes K 
corrected to redshift z = 0. We rest rict the sample to 15 < r < 
17.5, where r is the IPetrosianl d 19761) r-band extinction-corrected 
apparent magnitude. These bright and faint target limits define a 
complete sample (target selection is done only in r). The faint limit 
is r = 17.77 in some parts of the sky, but to be safe, we adopt 17.5. 
We also only consider galaxies with redshift 0.01 < z < 0.3. The 
low redshift cutoff excludes galaxies with large peculiar velocities. 
We primarily present results in the r band, but note that the g and 
i bands show very similar behavior. The u and z bands suffer from 
filter and chip problems, as well as low signal levels, and so we do 
not consider them. 

In order to isolate elliptical galaxies, we use the following 



three main criteria. Fi rst, the Sersic index in gri must satisfy n > 
2.5, in order to isolate lde Vaucoul eurs ( 1948) profiles (n = 4) over 
exponential profiles (n = 1). Second, we restrict the concentration 
index c = R90/R50 > 2.86 in r, where P50 and P90 are the Pet- 
rosian 50% and 90% light radii, respectively (for reasons explained 
shortly, we also require the concentration index c < 3.8 in gri). 
The criterion c > 2.86 se parates early and late-type galaxies with 
a com pleteness of ^82% dNakamura et al.ll2003l ; IShimasaku et al.l 
1200 lh . The concentration is not corrected for seeing effects, and 
thus a compact galaxy observed under bad seeing condition s will 
have an u nderestimated concentr ation (Bl anton et al.ll2003bl) . The 
sample of Nakamura et al. ( 2003) has eyeball classifications, how- 
ever, and the cutoff c = 2.86 was chosen to separate early and 
late types with minimal cross-contamination. For comparison, the 
Sersic index above does take into account the observed point 
spread function (PSF). In addition to the above cuts on the sur- 
face brightness profile, we also apply the colour criterion 01 (g- 
r) > 0.7, where 01 (g-r) is the g-r colour _K"-corrected to redshift 
O.l. lBlantonetalJ j2003b) found that this criterion isolates early- 
type galaxies from spirals. We note that the three criteria above 
are not completely independent. The most restrictive is the con- 
centration cut, followed by the Sersic cut and the colour cut. Their 
combination allows us to remove marginal outliers, especially SO 
galaxies, and thus focus on the normal elliptical population. 

In addition to the primary selection criteria above, several 
additional criteria are used to remove bad points, notably using 
pipeline warning flags, requiring a median signal-to-noise (S/N) 
> 10 in the spectrum, eliminating galaxies with very high or low 
(unphysical) values for a, requiring the de Vaucouleurs likelihood 
of the SDSS model fit to be greater than the exponential likeli- 
hood, and ensuring that all quantities have well-defined measure- 
ment errors. These quality restrictions remove roughly 7% of galax- 
ies identified by our elliptical criteria. For reasons explained in 
S |2.1| we also sepa rate out all BCGs identified in the C4 catalog of 
Miller et alj d2005t) from our main sample. The final sample (here- 
after known as the main normal elliptical galaxy sample) consists 
of 79,482 galaxies. 

The VAGC provides three different magnitude definitions for 
galaxies: a) a Petrosian magnitude, b) a magnitude derived from 
a Sersic fit, and c) a magnitude derived from a de Vaucouleurs fit. 
We elect to use Petrosian magnitudes, converted into total "Sersic"- 
like magnitudes usin g the concentration index c and the relations of 
iGraham et all d2005h for SDSS data. In particular, 



Mtot w M P -Piexp(c P2 ) 



(1) 



where Mtot is the total "Sersic"-like magnitude, Mp the Pet- 
rosian magnitude, and c the concentration. Pi = 4.2 x 10~ 4 and 
P2 = 1.514 for SDSS data. These transformations effectively re- 
cover the magnitudes that would be obtained from an ideal Sersic fit 
using only Petrosian quantities. The only assumption needed is that 
a galaxy's profile is well modeled by a Sersic fun c tion , although a 
Sersic fit is not actually used. The IGraham et alj d2005h relations 
also yield an effective radius R e (i.e., the half-light radius), using 



Re 



RSO 



1 - P 3 C P 4 ' 



(2) 



where R 5 q is the Petrosian half-light radius. P3 = 6.0 x 10~ 6 
and P4 = 8.92 for SDSS data. The above transformation equa- 
tions break down for c > 3.8, so we exclude such galaxies. This 
requirement removes less than 0.5% of galaxies, independent of 
luminosity, and does not effect our results. 

The conversion to an effective Sersic magnitude using Pet- 
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rosian quantities is desirable because (1) de Vaucouleurs mag- 
nitudes presume that all ellipticals are well-modeled as a 7? 1//4 
profile. In fact, however, the Sersic index ap pears to increase 
mildly with increa s ing luminosity (see, e.g., ICaon et al 1 1 19931 ; 
iGraham et ail 1 19961 ; iFerrarese et all [2006 ; we also see the same 
trend in our data; see i]3.2| l, and (2) Petrosian and Sersic quanti- 
ties themselves suffer from mild systematic errors. The fraction of 
total flux contained w ithin the Petrosian a perture is a function of 
the profile shape (see Graham et al. 2005 for an in-depth discus- 
sion). For example, the Petrosian aperture contains ~ 98% of the 
total flux for an exponential profile, but only ~ 80% for a de Vau- 
couleurs profile. This introduces a luminosity-dependent error into 
the total magnitude definition, given the mild increase in Sersic in- 
dex with increasing luminosity. Sersic magnitudes themselves do 
account for profile shape, and represent a total integrated luminos- 
ity. The re are, however, syste mati c problems with VAG C Sersic fits. 
As both lBlanton et al.|j2003) and lGraham et aljj2005h explain, the 
Sersic fluxes obtained from these fits are systematically underesti- 
mated, and become worse for increasing Sersic index n. At n = 5, 
the error is ~10%. Sersic fits are also more sensitive to the radial 
range over which the fit is carried out. 

For the reasons above, we elect to present results using the 
total Sersic-like magnitudes derived from Petrosian quantities. We 
note, however, that the entire analysis in this paper was repeated 
using the quantities obtained from the Sersic and de Vaucouleurs 
fits. We find results similar to those presented here, in particular, 
the same trends in the FP projections (see i]3.2t . 

The Sersic-like Petrosian quantities we use do not take into ac- 
count possible isophote ellipticity. The majority of galaxies in our 
sample, however, have axis ratios greater than ~ 0.6, and there is 
no strong dependence on abs olute mag nitude. As a check, we ap- 
plied a correction factor from West (2 005h to our adopted radii and 
repeated our analysis below. Our measured FP projection slopes 
changed by at most ~ 5%, with all the same trends being present. 
Ellipticity therefore does not have an important effect on our re- 
sults. 

In order to correct for passive evolution of early-type galaxies 
over the redshif t range probed by SDSS, we use the prescription of 
llin et al] d 19991) to correct all magnitudes; 

M r (z = 0) = M r {z) + Qz, (3) 

where M r (z — 0) is the r-band absolute magnitude corrected for 
evolution, M r (z) the observed absolute magnitude, z the redshift, 
and Q the evolution parameter i n magnitudes per redsh ift. We use 
the value of Q = 0.85 found bv lBernardi etail |2003b) for ellipti- 
cal galaxies in the r ba nd. Our net if +ev olution correction is very 
similar to that used by Wa ke et al] d2006h for luminous red galax- 
ies in SDSS in the r band, where K is the K correction at redshift 

2 = 0. 

All galaxies in the main sample have a Petrosian R50 > 1", 
larger than the typical seeing for the SDSS survey, which is a full 
width at half maximum of l'/4 in the r band. Thus none of the mea- 
sured radii are severely compromised by the PSF. As a quantitative 
check, we repeated our analysis of the FP projections imposing an 
additional cut that all galaxies have R50 > l"5, and found quanti- 
tatively similar results. 

We correct all velocit y dis persions for aperture effects follow- 
ing [jgrgensenetal] j 19951) and lwegner et al ] Jl999l) . This corrects 
all measured velocity dispersions to a standard circular aperture 
with a radius equal to one-eighth the effective radius of the galaxy: 



Ccorr — tineas ( ^ j , (4) 

where u mcas and a COTT are the measured and corrected velocity dis- 
persions respectively, i?flbcr the fiber radius from the SDSS (1.5), 
and R e the effective radius of the galaxy in arcseconds. We note 
that this correction assumes that all early-type gala xies have sim- 
ilar vel ocity dispersion profiles, regardless of Re. iBernardi et al] 
J2003dlbh provide a good discussion on this topic, but given that 
the correction depends very weakly on R e , these differences do not 
have any significant impact on our results. 

Finally, we define the effective surface brightness of a galaxy 
as follows, 7f -corrected and corrected for cosmological surface 
brightness dimming, 

^ = m + 2.51og 10 (27ri^) - 101og 10 (l + z) - K, (5) 

where m is the extinction- and evolution-corrected apparent mag- 
nitude, R e the effective radius (in arcseconds), z the redshift and 
K the Tf-correction at redshift z — 0. 

2.1 Brightest Cluster Galaxies 

As both lLauer et al] j2006tj) and IBernardi etail d2006h discuss, 
there is a problem with the standard SDSS photometry of lumi- 
nous BCGs (see also http://www.sdss.org/dr4/help/known.htmll. 
The problem is that the default sky subtraction removes the outer 
low surface brightness flux from the galaxy, resulting in an under- 
estimate of the luminosity and effective radii. We have indepen- 
dently come to the same conclusion. This problem is severe for 
BCGs because most of the luminosity is in low surface bright- 
ness emission. We were motivat ed to look into this p roblem by 
the lack of curvature in the initial lBernardi et alj d2003b) SDSS FP 
pr ojections at high luminosi ty, relative to the local BCG samples 
of Oegerle & Hoessel (1991) and Lauer (private communication & 
lLauer et al]|2006d) . 

We id e ntify BCGs in SDSS using the C4 catalog of 
iMiller et al] d2005t) . which is based on DR2 data, complete to 
z = 0.15. Clusters are defined in C4 on the basis of position, red- 
shift and colour, and the BCG is simply the brightest system in the 
cluster. In order to identify the C4 BCGs in the VAGC catalog, we 
use their declination and right ascension to find unique matches to 
within 1" in the DR4 catalog. We were able to successfully match 
744 out of 748 BCGs (some of which are duplicate objects), al- 
though 194 were too bright to be included in the spectroscopic sur- 
vey (i.e., brighter than the fiber saturation limit on SDSS). Of the 
remaining 550 BCGs, only 346 unique objects pass all our criteria 
in establishing an early-type galaxy sam ple, the most restricti ve be- 
ing the concentration cut. We note that Bernardi et al. I d2006h also 
find significant con tamination from non early-type galaxies in the 
IMiller et al] < |2005|) BCG catalog. The majority of this contamina- 
tion is likely due to a fiber collision restriction in SDSS data, such 
that the true BCG in a C4 cluster is missed in the s pectroscopic 
catalog dMiller et alj2005l ; lvon der Linden et al]|2006l) . 

To assess the accuracy of the standard pipeline photometry, 
we have reanalyzed images of a subset of high luminosity normal 
galaxies and BCGs. Sky subtraction was performed on each indi- 
vidual SDSS field containing a sample galaxy. All objects in the 
field were masked out and a t ilted plane was fi t to the remaining sky 
pixels. For more details see IWest et al] d2007l) . In the case where 
a sample galaxy crossed multiple fields, the fields were sky sub- 
tracted and then mosaicked. We then performed a 2-dimensional 
Levenberg-Marquardt least-squares de Vaucouleurs fit, masking 



^\ nnnn tj A c i\ yiMTi A c AAA nnft nf\n 



4 L. -B. Desroches et al. 



out nearby stars, galaxies, and the center to safely avoid the PSF 
core. 

Unlike the bulk of the normal elliptical galaxy population, we 
adopt a de Vaucouleurs fit instead of a Sersic fit for BCGs because 
of the subtle nature of BCG photometry. Sersic fits, with an extra 
degree of freedom, are sensitive to the range of the fit and tend 
to result in very large luminosities and radii for many BCGs with 
deep photometry (see £13 . 3b . Much of this flux, however, is probably 
better attributed to intracluster light rather than host galaxy light, 
though a proper decomposition into the two components is difficult 
and will likely require improved theoretical models. In the absence 
of a more appropriate choice, we adopt a de Vaucouleurs fit for the 
BCGs. 

Figure[T]shows a plot of the ratio of our fitted de Vaucouleurs 
effective radii to the catalog radii as a function of velocity dis- 
persion a for a sample of normal elliptical galaxies (black points) 
and BCGs (grey stars). We include a comparison to radii obtained 
from both the Sersic-like Petrosian measurements used in this pa- 
per and catalog de Vaucouleurs fits. For most of the normal galaxies 
the agreement is quite good, even for very massive and luminous 
galaxies with a ~ 300 km s -1 . This indicates that, for most nor- 
mal galaxies, the default catalog photometry is accurate. For the 
BCGs (gray stars), however, the catalog effective radii are smaller 
than our fitted values by factors of ~ 1.5-3. Similarly the total lu- 
minosities of these galaxies are underestimated in the VAGC by 
factors of sa 1.2-2. 5 . For t hese reasons, we exclude BCGs identi- 
fied by |Milleretai1 (2005) from our main sample, and reanalyze 
them separately. 

Despite the systematic underestimates in the catalog R e and 
luminosities for SDSS BCGs, BCGs represent only a small fraction 
of all elliptical galaxies and so their effect on our mai n sample as 
a who le is small. Indeed, only 346 BCGs identified by Mill er et all 
d2005h in DR2 are excluded from our main sample. As an indepen- 
dent check on this number, we note that using the number density of 
BCGs of 1.5 x 10' 5 h 3 Mpc~ 3 at z < 0.05 from lPostman & LaueJ 
we would expect w 400 BCGs in DR2, in reasonable agree- 
ment with the above numbers. For comparison, DR2 contains a to- 
tal of roughly 10000, 3100, and 400 elliptical galaxies that meet our 
criteria brighter than M r = -22.5, M r = -23, andM r = -23.5, 
respectively. Thus even if all BCGs were brighter than M r — —23, 
which they are not, BCGs would only become statistically signifi- 
cant (by number) for galaxies brighter than M r ~ —23. The pre- 
cise luminosity at which BCGs become statistically significant de- 
pends, however, on how the luminosity of a BCG is defined, i.e., 
how much of the extended low surface brightness emission is de- 
tected and how much is attributed to the galaxy rather than intra- 
cluster light (see i ]3.3t . 

To verify the small effects that BCGs have on our results for 
the normal elliptical galaxy sample, we performed our analysis of 
the FP projections on four data sets: a) DR4 with DR2 BCGs ex- 
cluded (our standard main sample); b) DR4 with no DR2 BCGs 
excluded; c) DR2 with DR2 BCGs excluded; and d) DR2 with no 
DR2 BCGs excluded. In all cases, our L(a) distributions and fits 
described below were virtually unchanged. The L(B e ) distributions 
and fits showed only small differences at the luminous end. These 
differences were of order 3%, roughly equal to our adopted la er- 
rors at the luminou s end. Because the identification of BCGs by 
iMiller et"al] d2005h is likely incomplete, we cannot rule out that 
contaminatio n by BCGs is presen t, particularly because we are us- 
ing DR4, and lMiller et al.l d2005h 's identification of BCGs is only 
available for DR2. Since DR4 contains roughly twice as many 
galaxies as DR2, we estimate that the number of remaining BCGs 



in our sample is similar to the number we have excluded (~ 350). 
The above comparisons between samples with and without DR2 
BCGs show, however, that although the photometry of some re- 
maining BCGs in our sample may be incorrect, they are not caus- 
ing significant errors in our analysis of the FP projections of the 
non-BCG elliptical population. 



3 FUNDAMENTAL PLANE PROJECTIONS 

3.1 Magnitude-Limited Sample Corrections 

Our sample is magnitude limited, as defined in Section [2] at both 
the bright and faint ends. In order to analyze the correlations be- 
tween observables for the elliptical population as a whole, we must 
correct for this selection effect. To accomplish this, we use the 
Vmax method, which weighs each point in the sample by 1/Knax, 
where V m &x is the maximum spatial volume one can observe a 
given galaxy in (adjusted for our redshift limits). We ignore K- 
correction and luminosity evolution in determining Vmax. They 
have nearly equal and opposite effects on the observed brightness 
of a galaxy (roughly one magnitud e per unit redshif t ), and thus do 
not aff ect our results significantly dShenetal.ll2003l ;i Blanton et al.l 
l2003al lcl). We choose the nonparametric V ma , x method over a max- 
imum likelihood analysis to avoid any assumptions about paramet- 
ric forms for the intrinsic distributions, even though we find them 
to be reasonably Gaussian. We note that the Vmax method is subject 
to inaccuracies at the faint end of the galaxy distribution, where an 
inhomogeneous and non-representative sample is given a very high 
weight. This well-kn own problem is ad dressed below. 

As discussed in ISheth et alj J2003h using a 1/V ma x correction 
represents the joint distribution of an observable X and luminosity 
L. The fit to this joint distribution is known as a bisector fit. In or- 
der to consider the distribution of the observable X at fixed L (i.e., 
< X | L >oc L), we must weigh each galaxy by l/[cf>(L)V max ], 
where cj>(L) is the luminosity function. The fit to this distribution 
is the inverse fit (i.e., X as a function of L). The inverse fit is more 
appropriate for these studies because the observable X is usually 
far more uncertain than the luminosity L. Treating L as an inde- 
pendent variable and assuming the scatter is largely due to X is 
thus the best course of action. 

One difficulty with the l/[</>(Z/)V m ax] correction is that a few 
galaxies are given high weight at both the faint and bright ends. We 
circumvent this issue by trimming our sample at the bright and faint 
ends. While this somewhat diminishes our results on the curvature 
in the FP projections, which are more pronounced at the extreme lu- 
minosity ends, it is necessary to prevent spurious fits dominated by 
small number statistics. Our criterion used to trim the sample is that 
we consider only galaxies with a weight Wi = l/[cj)(Li) V m ax,i] no 
larger than 10 times the mean, i.e., 



This effectively cuts out galaxies fainter than M ~ —19.9 and 
brighter than M ~ —23.7 in r band, removing a total of ~ 1.7% 
of galaxies. This cut is based on trial and error. Factor of few vari- 
ations about the choice in equation {6]l does not change any of our 
conclusions, while for a substantially more permissive cut the fits 
become overly dominated by a few galaxies, resulting in very un- 
physical fits. 
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3.2 Observed FP Projections 

Figure [2] shows the distribution of observables at fixed luminos- 
ity. The observables we consider are effective radius (R e ), velocity 
dispersion (<j), dynamical mass (Md yn oc R e <J 2 ), effective density 
(a 2 /R?), and effective surface brightness (p e ). We concentrate our 
discussion on the first two projections, but include all five for com- 
pleteness. In Figure[2] we separate the data points into small 0.25 
magnitude-wide bins and show the mean and la width of the dis- 
tribution in each magnitude bin; the distributions are well fit by 
Gaussians. 

Even by eye, Figure[2]shows pronounced curvature in the FP 
projections, with the exception of R e a 2 . We perform two tests to 
quantify this curvature. First, we perform a series of linear least- 
squares fits of observable X as a function of magnitude; we repeat 
these fits considering different subsamples with varying bright and 
faint end magnitude cutoffs. This quantifies how the slope of the 
FP projection depends on the luminosity range considered. Fig- 
ure [3] shows the resulting variation in the slope of the FP projec- 
tions, considering limits to the sample at the faint (black points) 
and bright (grey points) ends. The left-most black points and right- 
most grey points in any of the plots give the slope of the projec- 
tion over the entire sample of elliptical galaxies. The remaining 
points describe the slope of the FP projections over increasingly 
smaller sub-samples. The power-law slopes are defined as follows: 
R e oc L a , L oc a 13 , L oc (i? e cr 2 ) 7 , L oc (a 2 /R?) s , and fj, e oc U . 
The slope error is a statistical error using the measurement errors 
on a and R e . 

As an alternative way of quantifying curvature in the FP pro- 
jections, we consider a quadratic fit to the entire sample, of the 
form 



X = a + 6(log(L)) + c(log(L)) 2 , 



(7) 



where L is the luminosity and X is the observable. The parameters 
for these fits are presented in Table[T] In four of the five projections, 
the quadratic term is statistically non-zero. The observable R e a 2 is 
the only one without a clear detection of a quadratic term. Note that 
we are not ascribing any particular physical significance to this sec- 
ond order term, but merely using it to quantify curvature. By using 
the derivative of the quadratic fit with respect to log(L), we define 
the local slope of the L-X relation for all observables X. This is 
presented in Figure [4] The variations in these L-X slopes agree 
well with the trends seen in Figure[3] The statistical errors on these 
slopes are similar to those in Figure [3] for the largest subsamples 
and are omitted for clarity. 

To test the robustness of our local L-X slopes, we also per- 
formed 3rd order polynomial fits on all our relations, of the form 
X = a + b(log(L)) + c(log(L)) 2 + d(log(L)) 3 . The best fit local 
slopes of the L-X relations were nearly identical to our quadratic 
fits over the luminosity range of our sample, providing some confi- 
dence in the robustness of the results in Figure|4] In the case of the 
FJ relation, however, a quadratic expansion with respect to log(L) 
might not be the most appropriate description of the data, given the 
observed sign of the curvature in Figures [2] & [3] In equation 0, 
a — > as L — » 00 (because c < in Table[TJ, which is clearly un- 
physical. An a priori more reasonable expansion of the FJ relation 
would be log (L) = a + bX + cX 2 , where X — log(cr) (see also 
IWvithell2006l who carried out a similar fit for Mbh-o"). This now 
has the correct asymptotic behavior at large a. The resulting best 
fit local L-a slopes are again almost identical to the results shown 
in Figure|4]over the range in L we are interested in. 

It is clear from Figures [3] and [4] that substantial curvature ex- 



ists in the FP projections across the entire sample of normal (non- 
BCG) elliptical galaxies, with the exception of the dynamical mass 
vs. luminosity relation which is well-fit by a single power-law with 
L oc MJ , with 7 « 0.86, which implies M dyn /L oc L 016 . 
This d ynamical mass-to-light ratio agrees well with the DR 2 re- 
sults o flPadmanabhan et al.U2004l) . |Padmanabhan et alj fc004h also 
show that the stellar mass-to-light ratio M+/L of early type galax- 
ies is constant with lum inosity, using stellar masses determined by 
iKauffmann et alj d2003h . Thus 7/I accounts for at least some of 
the tilt of the FP in the r band relative to the virial theorem expec- 
tation. An increasing dark matter fraction with increasing stellar 
mass coul d be responsible for this d ynamical mass-to-light ratio 
variation (Bovlan-Kolchin e t al.ll2005T) . 

As a check on the importance of non-homology for tilting the 
FP, Figure [5] shows the Sersic index n as a function of luminosity 
for our sample. Although the range of n at any given luminosity 
is quite large, there is a clear trend of increasing n for more lumi- 
nous galaxies, as has been found by other authors iCaonetalJl993l : 
iGraham et al.ll996l ; lFerrarese et al.l2006h . This trend is slightly un- 
derestimated in the VAGC due to the systematic problems associ- 
ated with the Sersic fits, which result in fluxes, ra dii, and Sersic 
indices being underestimated by ~10% at n = 5 ( Blant orTet al.1 
20051). The variati on in n we find is significantly less than what 



CiottietalJ jl996l) find is required for non-homology to fully ac- 



count for the tilt in the FP across the entire luminosity range of 
early type galaxies. It is therefore likely that non-homology, while 
present, does not fully account for the tilt in the FP. 

Aside from L vs. Md yn , all of the other FP projections ex- 
hibit a local slope that varies systematically with luminosity from 
the faint end of the sample to the bright end. The slope a of the 
radius-luminosity relation for elliptical galaxies varies systemati- 
cally with luminosity, from a ~ 0.5 at M T w —20 to a w 0.7 
at M r m —24 ; for the en t ire sa mple , we find a fa 0.6, in good 
agreement with lShen etall d2003h and lBernardi et alj d2003bh . Our 
slope at faint lumin osity of a ~ 0.5 agrees with that found by 
lLauer et al] d2006al) for their lower-lumino sity ellipticals. At the 
luminous end, however, lLauer et alj j2006j) find a ~ 1.1, much 
steeper than we se e, though their samp le consists primarily of 
BCGs in the Vband. lBernardi et alj ( ]2006f) attribute curvature in the 
radius-luminosity relationship at high luminosity in the full SDSS 
sample to an increasing fraction of BCGs dominating the sample. 
Our results show, however, that the normal non-BCG population 
exhibits a steepening at high luminosity that does not appear to be 
attributable to contamination by BCGs (see <i2. 1 b . 

The slope /3 of the FJ relation also steepens from ~ 3 
at M r « -20 to » 4.5 at M r ~ -24; the canonical 
FJ slope of L oc a 4 is found in the middle of our sample, at 
M r = —22.5. Our slo pe at fa i nt lum inosity for the FJ relation is 
si milar to that found by Tonrvl (fl98 lh . but not as flat as that found 



bv lMafkovic & Guzmar J2005h . who report L 



. a .ui±u. 30 for a 

sample of faint early- type galaxies; their sample, however, extends 
to much fainter magnitudes (—22.0 < M r < —17.5 mag). At the 
luminous end, our value of fj is si gnificantly smaller than the value 
of (3 w 8 found bv lLauer et afl2 006a. although this is again a result 
of their sample being primarily BCGs. 

The density-luminosity and surface brightness-luminosity re- 
lations shown in Figures |2j [4] indicate that low-luminosity ellip- 
ticals have roughly constant densities and surface brightnesses, 
while more luminous ellipticals become progressively less dense 
and have lower surface brightness with increasing luminosity. 

As discussed in Sj2] we find similar results for the FP pro- 
jections regardless of the magnitude we use to define the sample 
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(Petrosian, Sersic, or de Vaucouleurs). More specifically, the lo- 
cal slopes of the quadratic fit to the FJ relation are nearly identical 
for all three magnitude types. For the radius-luminosity relation, de 
Vaucouleurs fits result in a ~ 0.45 at the faint end and a ~ 0.75 
at the bright end, whereas with Sersic fits the slope ranges from 
a ~ 0.5 to a ~ 0.8, both similar to our default analysis (Fig.[4]l. 

3.3 Brightest Cluster Galaxies 

BCGs are at the extreme end of the luminous elliptical galaxy popu- 
lation; it is thus interesting to compare our results on the FP projec- 
tions for luminous ellipticals to those of BCGs. Figure[6]shows our 
non-BCG SDSS data for the Re-L correlat ion of elliptical s , alon g 
with a local sample of BCGs studied by Gonz alez et al. I J2005t) . 
who extend previous BCG studies by studying the very extended 
low surface brightness emission at large radii and comparing a vari- 
ety of full two dimensional models for BCG s urface brightness pro - 
files. We also include the local BCG sample o flLauer et aljj20 06b). 
We use SDSS / band data shifted to the I c band of lGonzalez et al.l 
using the median i-z colour and the transformation equations ob- 
tained from the SDSS website. (For ease of comparison, we also 
include the approximate r band magnitudes on the upper axis, us- 
ing the med i an colo ur of r-I c = 0.87 for the SDSS sample.) The 
lLauer et al] d2006bl) points are shown as dots, shifted to the I c 
band using a median colour of V -I = 1.4 for this sample. The 
solid grey points in Figure [6] are iGonzalez et alt s 1-component 
de Vaucouleurs fits for R e and L. They also argue, however, that 
Sersic fits and fits with two de Vaucouleurs components (an in- 
ner and outer component) provide a much better description of the 
surface brightness profiles of BCGs; the resulting R e and L are 
shown by stars and open circles in Figure [6] respecti vely. F or the 
two-component de Vaucouleurs fits. IGonzalez et al. (2005) argue 
that the inner component may correspond to a "normal" elliptical 
galaxy while the outer component represents an extended envelope. 
Figure|6]shows, however, that most of the effective radii for Gonza- 
lez's inner component fits lie well outside the observed distribution 
of R e for normal ellipticals. This, together with the large scatter in 
the Re-L relation for these inner component fits, makes it some- 
what difficult to physically interpret the two component model for 
BCG surface brightness profiles. 

The 1-componen t de Vaucouleurs fits for R e and L from 
IGonzalez et alj d2005h yie ld an Re-L correlation of R e oc 
L 1±01 . | Lauer et al] J 2006a) find a similar result for BCGs, as do 
iBernardi et al] ( 120061) . For comparison, we find that R e oc L ' 8 at 
the luminous end of the normal elliptical galaxy population in the 
i band (note that this is slightly steeper than our r band slope). In 
Figure [7] we show the Faber- Jac kson relation with ou r SDSS sam- 
ple and the local BCG sample of Lauer et al. (2006a), both shifted 
to the I c band as described above. Lau er et al. I d2006al) find L oc a 
for local BCGs, while we find L oc cr 4 ' 5 for the most luminous el- 
lipticals. The Faber- Jackson and R e -L scalings for BCGs are thus 
significantly steeper than those for normal elliptical galaxies, even 
at the luminous end of the latter population. This is consistent with 
significant growth of BCGs via dissipationless mergers (see Sj4] also 
iBovlan-Kolchin et al]|2006l) . 

In Figures [8] & [9] we compare the /--band radius-luminosity 
and FJ relations of our non-BCG main sample to our reanalysis 
of the C4 BCGs. The BCGs were r eanalyzed as desc ribed in §2. 11 
using the sky subtraction method of lWest et al] d2007l) . The radius- 
luminosity relation for BCGs appears distinct from the non-BCGs, 
exhibiting a steeper slope of a ~ 1.1. By contrast, the FJ rela- 
tion for BCGs appears to smoothly match the non-BCG sample, 



although it flattens significantly above M r = —23. The overall 
FJ slope for the BCGs is f3 ~ 4.2, similar to the slope at the lu- 
minous end of the non-BCG sample. Brighter than M r = —23, 
however, the FJ slope for BCGs is [3 ~ 5.9, while fainter than 
M r = —23, the slope is /3 « 3.0. These results agree rea- 



sonably well with Bernardi et al. (2006) and lLauer et"alT (l2006a), 
but Ivon der Linden et all J2006h find a much shallower radius- 
luminosity relation of R e oc L°' 65±0 02 for BCGs, consistent with 
normal elliptical galaxies. The latter sample is based on a careful 
reanalysis of C4 clusters, but with definitions of radius and lumi- 
nosity that differ from those used here. 

Most of the luminosi ty in BCGs arises in e xtended low surface 
brightness emission (e.g jGonzalez et alj2005h . Thus observations 
of different depths will yield different luminosities and the inferred 
"galaxy" luminosity and radius hinge critically on how much of 
the extended low surface brightness emission is attributed to intra- 
cluster light rather than gala xy light. The importance of this is il- 
lustrated by IGonzalez et alj 's Sersic fits in Figure [6] which give 
an Re-L correlation of R e oc L 1 - 8 "' 2 , very different from the 
R e oc L 1 correlation of their one component de Vaucouleurs 
fits. This highlights the subtlety in defining the radii and luminos- 
ity of BCGs, which needs to be taken into account when comparing 
the properties of normal ellipticals and BCGs. 



4 DISCUSSION 

We have analyzed the fundamental plane projections of elliptical 
galaxies as a function of luminosity, using a sample of ~ 80, 000 
galaxies drawn from DR4 of the SDSS. We have separated BCGs 
from our main elliptical galaxy sample and reanalyzed their pho- 
tometry. The observables we consider are effective radius (Re), ve- 
locity dispersion (cr), dynamical mass (Afd yn oc R e a 2 ), effective 
density (a 2 /R^), and effective surface brightness (fi e ). With the ex- 
ception of the L-Mdyn correlation, we find clear evidence of varia- 
tions in the slope of the FP projections with luminosity (Figs.[3]&[4]l 
in the normal elliptical galaxy sample. The trends we find are that 
the radius-luminosity and Faber-Jackson relations are steeper at 
high luminosity relative to low luminosity and that more luminous 
ellipticals become progressively less dense and have lower sur- 
face brightnesses than lower luminosity ellipticals. These trends are 
consistent with existing results in the literature, though our results 
have much better statistics. It is interesting to note that the results 
we find appear to continue to even lower luminosity spheroidal sys- 
tems. In particular, for a sample of dwarf elliptical (dE) and dwarf 
spheroidal (dSph) galaxies in the B band, the radius-luminosity 
power-law slope a sa 0.28-0.55 and the Faber-Jackson power-law 
slope (5 « 1.5-2.5 dde Riicke et"aHl2005h . broadly consistent with 
an extrapolation of our results to eve n lower luminosities . Thes e 
trends are also seen in the dE sample o flGraham & Guzmanl J2003h . 
We note, however, that there are likely structural differences be- 
tween luminous ellipticals and dEs that complicate this simple ex- 
trapolation. 

The variations in the FP projections we find at the luminous 
end are consistent with less and less dissipation during the forma- 
tion of elliptical galaxies with increasing luminosity. Less dissi- 
pation would result in increasingly larger, less dense, lower sur- 
face brightness, and lower a galaxies, relative to fiducial power- 
law scalings for the elliptical galaxy population. This interpreta- 
tion is consistent with a variety of other observational evidence 
for less dissipation in the formation of luminous ellipticals (e.g., 
iKormendv & Bender|[l996l) . 
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Luminosity (or galaxy-mass) dependent variations in the FP 
projections should in principle provide a strong constraint on 
galaxy formation models. For example, numerical simulations of 
gas-free ("dry") merger remnants find that the FP projections have 
a ~ 0.7 and /3 ~ 4 for very wide orbits, while lower angular 
momentum and higher energy orbits have less energy exchange be- 
tween stars and dark matter and thus the FP proj ections steepen to 
a ~ 1 and (3 ~ 6 CBovlan-Kol chin et alj|200q) . These values of 
a and ft are consistent with the trends we find in the data at the 
luminous end, suggesting an increasing importance of dry mergers 
for luminous ellipticals. For lower-luminosity ellipticals, however, 
dry mergers are less likely to be important in determining the prop- 
erties of elliptical galaxies. Current simulations of elliptical galaxy 
formation via mergers of gas-rich disks find that the remnants lie 
roughly on the FP and its projections, with the gas fraction of the 
progenitor galaxies strongly affecting the propertie s of the resulting 
spheroidal merger remnants (Roberts on et al.l2006h . These simula- 
tions also find more dissipation in the formation of lower lumi- 
nosity ellipticals, which is consistent with our data. It remains to 
be seen whether simulations of gas-rich mergers can also account 
for the systematic variations in the FP projections with luminosity 
(e.g., with a varying gas fraction). 

We have also reanal yzed the pho t ometr y of SDSS BCGs using 
the C4 cluster sample of lMiller et al.l |2005); the standard pipeline 
photometry underestimates the luminosities and radii for luminous 
BCGs (Figure Q3 due to an overestimate of sky background. With 
our corrected photometry, we find that the radius-luminosity rela- 
tion of BCGs is noticeably steeper than that of normal (non-BCG) 
ellipticals and the FJ relation for BCGs flattens at the luminous 
end. The same tren ds are also seen in local samples of BCGs (eg., 
Laueretal . 2006a). These observational results are consistent with 
significant growth of BCGs by dissipationless mergers. 
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4.1 Massive Black Holes 



For the bulk of the elliptical galaxy population, both the Mbh-o" 
and Fabe r- Jackson rela tions scale roughly as a 4 . As noted in §1, 
however, Wyithe (2006) argued that the A/bh-o" relation steepens 
at high black hole masses and flat tens at low masses, albeit with 
only modest statistical sig nificance. lGreene & Hold2006h see a sim- 
ilar flattening at low black hole masses and lLauer et alj d2006al) ar- 
gue for a similar steepening at high black hole masses based on 
the properties of the surface brightness cores in luminous ellipti- 
cal galaxies. These luminosity-dependent variations in the Mbh-0" 
relation are qualitatively similar to the luminosity dependent vari- 
ation in the FJ relation. Together with a constant M+/L, these re- 
sults imply a roughly constant v alue of Mbh/M+ ac ross a range 
of early type galaxies (see, e.g., lHaring&Rixll2004l for a direct 
compilation of Mbh/M* that bears this out). For massive galaxies, 
the approximate constancy of Mbh/M* may be a consequence of 
dissipationless mergers, which preserve Mbh/M+ in the absence 
of ene rgy loss by gravitational w aves during black hole coales- 
cence ( Bovlan-Kol chin et alj[2006r) . For lower luminosity ellipti- 
cals, however, gas dynamics is important and a constant value of 
Mbh/M* implies a constant relative efficiency for forming stars 
and massive black holes in a given gravitational potential well. 
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Table 1. Quadratic Fitting Parameters 



FP Projection 


a 


b 


c 


log(_R e ) vs. log(L) 


1.50 ±0.19 


-0.802 ± 0.044 


0.0805 ± 0.0026 


log(o-) vs. log(L) 


-1.79 ±0.23 


0.674 ± 0.053 


-0.0234 ± 0.0030 


log(7J e <7 2 ) vs. log(L) 


-3.37 ± 0.59 


0.81 ±0.14 


0.0199 ± 0.0078 


log((T/R e ) 2 vs. log(L) 


-4.82 ± 0.75 


2.54 ±0.17 


-0.1845 ±0.0098 


lie VS. log(L) 


43.4 ± 2.2 


-5.77 ±0.51 


0.350 ± 0.030 



All fits are X = a ± b(log(L)) + c(log(L)) 2 , where log(L) is the lu- 
minosity (M = — 2.51og(L)) and X is the observable. L is calibrated to 
the AB system, in which a magnitude object has the same counts as a 
F v = 3631 Jy source; a is in units of km s~ 1 ; R e is in units of kpc; fi e is 
in units of mag arcsec~ 2 . 



West, A. A., et al. 2007, AJ, submitted 

Wyithe, J. S. B. 2006, MNRAS, 365, 1082, Erratum: Wyithe, 
J. S. B. 2006, MNRAS, 371, 1536. 
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Figure 1. A comparison of our fitted R e to the catalog R e as a function 
of <r. R e fl t is the effective radius fro m a de Vaucouleurs fit to a reanalysis 
of SDSS photometry using West et ail 's (2007) sky subtraction. i?, e , pc tro 
is the effective radius derived from catalog Sersic-like Petrosian measure- 
ments, as described in J5] R e <j c v i s me effective radius derived from cat- 
alog de Vaucouleurs fits. Black points are a random sub-sample of high-L 
normal ellipticals. Grey stars are a random sub-sample of high-L BCGs. 
The BCG effective radii (and luminosities) are consistently under-estimated 
with standard SDSS photometry, while most normal elliptical galaxies ap- 
pear fine. Two BCGs have R e fit/Re,pctro = 4.0; these are omitted for 
clarity. 
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Figure 2. Distributions of effective radius R e , velocity dispersion a, dy- 
namical mass R e a 2 , effective density (a/R e ) 2 , and effective surface 
brightness fj, e in 0.25 mag wide bins, for a sample of 79,482 non-BCG 
early-type galaxies drawn from SDSS DR4 that pass concentration, Sersic, 
and colour cuts. Points denote the peak of a gaussian fit to the distribution 
in a given bin, while the error bars correspond to the la width of the distri- 
bution. 
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Figure 3. Slope of power-law fits for five fundamental plane projections 
of our non-BCG elliptical galaxy sample as a function of r-band magni- 
tude cutoff. Black points represent samples with a lower magnitude cutoff 
-^cutoff, where the sample is increasingly restricted to high luminosities. 
Grey points represent samples with an upper magnitude cutoff M cuto fj, 
with the sample increasingly restricted to low luminosities. These fits show 
that, with the exception of L vs. dynamical mass (R e a 2 ), the slopes of the 
FP projections depend on the luminosity of the elliptical galaxy sample. 



nnnn D A C A yIMD A C AAA AAA AAA 



12 L.-B. Desroches et al. 





-20 -22 -24 

M, [mag] 



-20 -22 -24 

M, [mag] 





-20 



-22 
M, [mag] 



-24 



Figure 4. Local slope of the fundamental plane projections of our non- 
BCG elliptical galaxy sample as a function of r-band magnitude. The data 
in Figure[2]are fit to a quadratic function of the form X = a + f>(log(L)) + 
e(log(L)) 2 and the local slope is defined by the derivative of the fit with 
respect to magnitude. Fit parameters are given in Table[T] These results for 
the variations in the fundamental plane projections with luminosity are very 
similar to those of Figure [3] The statistical errors on the slopes are similar 
to those in Figure[5]for the largest subsamples and are omitted for clarity. 
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Figure 5. Distributions of Sersic index in 0.25 mag wide bins for our non- 
BCG elliptical galaxy sample. Points denote the mean of a gaussian fit to the 
distribution in a given bin, while the error bars correspond to the lcr error on 
the mean. The mild variation in index with luminosity for low-luminosity 
ellipticals implies a non-homology within the elliptical galaxy population. 
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Figure 6. Comp arison between our non -BC G SDSS sample and the local 
BCG samples of iGonzalez et al l J2005l) and lLauer et"afl J2006bt) in the J" c 
band. The corresponding r band magnitude is estimated using the median 
colour of r-I c = 0.87 for the SDSS sample. SDSS data (black circles) are 
the effective radii vs. luminosity distributions in the ;' band, with BCGs ex- 
cluded, shifted to the I c band of Gonzalez using a median i-z colour of 0.22 
and the transformation equations obtained from the SDSS website. Also 
shown are Gonzalez's one component de Vaucouleurs fits for the effective 
radius and luminosity of local BCGs (grey circles), the inner component of 
Gonzalez's two component de Vaucouleurs fits to the photometry of BCGs 
(open circles), and the radii and luminosities from Sersic fits (stars). The 
lLauer et alj <2006bt) BCGs (dots) are one component de Vaucouleurs fits, 
and have been shifted to the I c band u sing a median colour o f V-I = 1.4. 
The two de Vaucouleurs samples from Gonzalez et al. 12005) can be fit by 
R e oc ^1-0±°'1 > which is steeper than the R e tx L o s correlation of nor- 
mal elliptical galaxies at the bright end in band. The Sersic sample is even 
steeper, with R e oc L 1 - 8 ^ - 2 . 
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Figure 7. Com parison between our non-BCG SDSS sample and the local 
BCG sample of lLauer et all (2006a) in the I c band. The corresponding r 
band magnitude is estimated using the median colour of r-I c = 0.87 for 
the SDSS sample. SDSS data (black circles) are the velocity dispersion vs. 
luminosity distributions in the i band, with BCGs excluded, shifted to the 
I c band using a median i-z colour of 0.22 and the transformation equations 
obtaine d from the SDSS website. The BCG data points from lLauer et al J 
l2006ah have been shifted to the I c band using a median colour of V-I = 
1.4. 
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Figure 8. Comparison of R e vs. L for our non-BCG SDSS sample and the 
C4 BCGs, whose photometr y has been reanalyzed using the sky subtraction 
method of IWest et al] <2007t) . The BCGs are disjoint from the non-BCGs, 
with a larger R e -L slope, and do not smoothly transition into the non-BCG 
relation. 



nnnn n a c a tmtj a c aaa fwi nnn 



Luminosity Dependence in the Fundamental Plane Projections of Elliptical Galaxies 



T 



2.6 - 



2.4 

in 

e 



° 2.2 



2.0 - 



• SDSS non-BCG 

• SDSS BCG 



-19 



-20 



-21 



-22 
M f [mag] 



-23 



-24 



-2 



Figure 9. Comparison of a vs. L for our non-BCG SDSS sample and the 
C4 BCGs, whose photometry has been reanalyzed using the sky subtraction 
method of IWest et alj J2007I) . The BCGs do not clearly differentiate them- 
selves from the non-BCGs, although they exhibit strong variations in slope 
at the high luminosity end. 
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